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Reflection of light from an anisotropic medium 



Abstract 



We present here a general approach to treat reflection and refraction of light of ar- 
bitrary polarization from single axis anisotropic plates. We show that reflection from 
interface inside the anisotropic medium is accompanied by beam splitting and can create 
surface waves. 

1 Introduction 

Theory of anisotropic media in optics alienates by its complexity, phenomenology and not 
transparent physics [1-5]. The most recent paper related to this field is published in Am. J.Phys.- 
[6] only in 1977. It is understandable that no one likes to discuss this topic. In the science 
of almost two hundred years old it seems hopeless to propose something new which can be 
attractive to students. Nevertheless, in this paper we will do just that. We will present an 
approach similar to that, used for description of elastic waves in anisotropic media [7] 

The single axis anisotropic nonmagnetic media in electromagnetism can be described by a 
matrix of dielectric permittivity e with matrix elements [1] 



where eo is isotropical part, and anisotropy is characterized by the unit vector a with compo- 
nents di and by anisotropy parameter e'. With this matrix we will first consider propagation of 
electromagnetic waves in a homogeneous anisotropic medium, then refraction at the interface 
between anisotropic and isotropic media, and after that discuss reflection and transmission of 
an anisotropic plain plate of finite thickness. Time to time along the text we make digression 
to isotropic media, because on one side it helps to check correctness of our formulas, and on 
the other side it is useful, because isotropic media also represent some difficulties. 

Of course, some our formulas look huge, but we don't care about it. We do not even derive 
some formulas and leave it to those who will need it. In fact, it is sufficient to know how to 
derive them, and if necessary to use a computer. 

2 Waves in an anisotropic medium 

The wave equation for, say, electric field E(r,t), is obtained from Maxwell equations. In a 
homogeneous nonmagnetic (/i = 1) anisotropic medium it is 



eo$ij + e'didj, 



(1) 




(2) 



Solution of this equation can be accepted in the form of a plain wave 



E(r,t) = £exp(ik ■ r — iut), 



(3) 



where £ is a polarization vector, which can be of not unit length. 
After substitution of (3) into (2) we obtain 

k 2 £ - k(k ■ S) = k 2 e£, 



(4) 



where ko = u/c, and e£, according to (1) is 

e£ = e £ + e'a(a ■ £). (5) 

In isotropic media polarization vector £ can have arbitrary direction perpendicular to the 
wave vector k. In anisotropic media it is not so. 

Besides (5) the field E(r,t) should also satisfy Maxwell equation 

V • eE(r, t) = 0, (6) 

from which, after substitution of (3), it follows a limitation on possible choices of directions for 
£: 

e (k-£) + e'(k-a)(a-£) = 0. (7) 

If k is not parallel to a, we have three independent vectors a, n = k/k and e± = [a x n], 
which can serve a basis of a coordinate system. The basis is not orthonormal, nevertheless the 
vector £ can be represented in this basis as 

£ = aa + /?K + 7ei (8) 

with some coordinates a, (5 and 7, which are not arbitrary, as will be now shown. 
Substitution of (8) into (7) gives 

e [k(3 + (k ■ a)a] + e'(k ■ a)[a + P(k ■ a)] = 0, (9) 

from which it follows that 

* (K ' a)(1 + ^a, (10) 



1 + t](k ■ a) 2 

where 77 = e'/e . Substitution of (10) into (8) gives 

/ (k-o)(1 + t])\ , . 

£ = a [a-K \ J ; ^ + 7 ei = ae 2 + 7ei, H 

which shows that in fact we have only two independent vectors for expansion of £: vector 
e± — [a x k] orthogonal to the plane of a, k, and the vector 

e 2 = a-n— — —. (12) 

1 + t](k ■ a) 2 

It is worth to note that e 2 is not an eigen vector of matrix e, but the matrix e transforms 
e 2 into a vector orthogonal to k 

ee 2 = €i(9)[k x [a x «]], (13) 

where 

e n°) = e o— — - ( ^ = e o^ Tfl, (14) 

l + r](K-a) 2 1 + r] cos 2 9 

and 9 is the angle between vectors a and k. It is seen that, if the wave propagates along a, 
the dielectric permittivity becomes ei(0) = eo, and, if the wave propagates perpendicularly to 
a, the dielectric permittivity becomes ei(7r/2) = eo + e'. 
With account of (14) we can represent (12) as 

e 2 = a - k(k ■ a)ei(9)/e . (15) 
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In the limit 77 — >■ 0, when the medium becomes isotropic we obtain 

lime 2 = [re x [a x re]]. (16) 

Now we will show that in general a plain wave in anisotropic media can have polarizations 
only along e 2 or e.\. 

Indeed, substitution of (11) into (4) with account of (13) gives 

k 2 [ae 2 + 7ei - are(re • e 2 )\ = k 2 {ae x (6) [re x [a x re]] + 7e ei). (17) 
Multiplying it with e\ we obtain 

{k 2 - fc 2 e ) 7 e? = 0. (18) 
Therefore, if 7 7^ 0, (18) can be satisfied only, when 

k 2 = k 2 e . (19) 

Multiplying (17) with e 2 we obtain 

(k 2 -k 2 e 1 (6))a[l-(a-K) 2 } = 0. (20) 
Therefore, if a 7^ 0, and a 7^ re, (20) can be satisfied only, when 

fc 2 = fcgexP), (21) 

where ei(0) is given in (14). Since the length of fc is different for two polarization vectors, 
therefore a single plain wave can exist only with a single polarization along either e 2 or e.\. 

We will call "transverse" the mode with polarization £1 — e±, and "mixed" the mode with 
polarization along £ 2 = e 2 . The mixed mode contains a longitudinal component - a component 
along the wave vector k. We think that such a nomenclature is better than common names: 
"ordinary" for S\ = ei, and "extraordinary" for £ 2 = e 2 , because our names point to physical 
peculiarities of these waves. 

2.1 Magnetic fields 

Every electromagnetic wave besides electric contains also magnetic field. For simplicity we 
assume that n — 1. From the equation V ■ H = 0, which is equivalent to k ■ H — 0, it follows 
that the field H is orthogonal to k. It is also orthogonal to £, which follows from the Maxwell 
equation 

-\yxE(r,t)] = -^H(r,t). (22) 

After substitution of (3) and 

H(r,t) = Hex.p(ik-r -iut), (23) 
with polarization vector H of the field H, we obtain 

H = ±[Kxe\, (24) 
k 

where k Q = u/c, and re = k/k is a unit vector along the wave vector k. For transverse and 
mixed modes, respectively, we therefore obtain 

Hi = —[re x ei] = —[re x [a x re]], "H 2 = —[re x e 2 ] = —[re x a], (25) 
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and the total plain wave field looks 

^/(r, t) = ipj exp(ikj ■ r — iut), (26) 

where ipj = Sj + "Hj, and j denotes mode 1 or 2. In isotropic media we also can choose, say 
£ = [a x k] and H — [k x [a x k]\. However, there a can have arbitrary direction, therefore 
the couple of orthogonal vectors £ and H can be rotated any angle around the wave vector k. 



3 Reflection from an interface with an isotropic medium 

Imagine that our space is split into two half spaces. The part at z < is an anisotropic medium, 
and the part at z > is vacuum with 6q = 1, rj = 0. We have two different wave equations 
in these parts, and waves go from reign of one equation into reign of another one through the 
interface where they must obey boundary conditions imposed by Maxwell equations. 

Let's look for reflection of the two possible modes incident onto the interface within the 
anisotropic medium. 



3.1 Nonspecularity and mode transformation at the interface 

First we note that reflection of mixed mode is not specular. Indeed, since direction of k after 
reflection changes, therefore the angle 9 between a and k does also change, and k, according to 
(21), changes too. However the component k\\ parallel to the interface does not change, so the 
change of k means the change of the normal component k±, and this leads to nonspecularity of 
the reflection. 

Let's calculate the change of k± for the incident mixed mode with wave vector k 2r , where 
the index r means that the mode 2 propagates to the right toward the interface. For a given 
angle 9 between k 2r and a we can write 



k 



2r_L 



e k 2 (l + r]) 

\ 1 + Tj COS 2 9 



- kL (27) 



however the value of k 2r ± enters implicitly into cos6>, so to find explicit dependence of k 2r ± on 
a it is necessary to solve the equation 

k\ + x 2 + rjik^l ■ a) + x{n ■ a)) 2 = fcjje (l + rj), (28) 

where x denotes k 2r ±, n is a unit vector of normal, directed toward isotropic medium, and I is 
a unit vector along fey, which together with n constitutes the plane of incidence. Solution of 
this equation is 



-77fc||(n • a) (I ■ a) + ^e k 2 (l + 77) (1 + r)(n ■ af) - Af (1 + r](l ■ a) 2 + V (n ■ a) 2 ) 

X = 1 l+^^-a) 2 ^ • (29) 

The sig n chosen b efore square root provides the correct asymptotics at i] = equal to isotropic 
value y^o&o — k 2 . 

In general vector a is representable as a = an + (31 + 7*, where t = [nl] is a unit vector 
perpendicular to the plane of incidence. The normal component k 2r ± depends only on part of 
this vector a' = an + (31, which lies in the incidence plane. If we denote a = \a'\ cos(# a ), (3 = 
I a' I sin(# a ), where \a'\ is projection of a on the incidence plane, and introduce new parameter 
rj' = i]\a'\ 2 < 77, then formula (29) is simplified to 

-rfk\\ sin(2fl a ) + 2^0^(1 + 77) [1 + rf cos 2 (# a )] - k 2 {\ + rf) 
hr± ~ 2[l + 7/cos 2 (^)] • (30) 
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Figure 1: Variation of change A of normal components for reflected and incident waves in 
dependence of angle 9 = 9 a of anisotropy vector with respect to normal n. Vector a is supposed 
to lie completely in the incidence plane. The curves Aij represent dimensionless ratio Ajj(6>) 
given by (34), (35) and (36). The curves were calculated for rj — rf' — 0.4 and q = k\\/koy/eo = 
0.7. 




For the reflected mixed mode (mode 2, propagating to the left from the interface) an equation 
similar to (28) looks 

fcj + x 2 + 7]{k\\ {I -a)- x(n ■ a)) 2 = fcge (l + rf), (31) 
where x = k 2 i±, and its solution is 

rfk\\ sin(2fl a ) + 2^0^(1 + 77) [1 + rf cos 2 (# a )] - k 2 (l + rf) 



"211. 



2[l + Vcos 2 (^)] 



(32) 



We see that the difference of the normal components of reflected and incident waves of mixed 
modes k 2 i± — k 2r ± is 

rfkn sm(29 a ) 

k 2 i± - k 2r± = — — — — , (33) 

1 + cos 2 (9 a ) 

In the following we will present such differences in dimensionless variables 

k 2i ± - k 2r± _ v'qsH^a) + 2y/(l + r ? )[l + Vcos 2 (^ a )]-g 2 (l + if) 



±22 



k ^~o 2[l + Vcos 2 (^ a )] 



(34) 



where q 2 = k 2 /k 2 e . The reflection angle depends on orientation of anisotropy vector a and it 
can be both larger than the specular one, when 9 a > 0, or smaller, when 9 a < 0. 

In the case of transverse incident mode the length k — \k\ of the wave vector, according to 
(19), does not depend on orientation of a, therefore this wave is reflected specularly. 

Every incident mode after reflection creates another one, without another mode it is im- 
possible to satisfy the boundary conditions. Let's look what will be the normal component of 
the other mode. If the incident is of mode 2, reflected transverse mode (mode 1 propagating 



to the left, away from the interface) will have ku± = y^o^o — k 2 . Therefore according to (30) 
the difference A i2 = (ku± — k 2r ±)/k ^/eo is 



-))'<)sm(2« a )+2j(l + ,,)[l+r,'co S 2 (ej]- 9 2 (l + ,,') 

2|1 + ,W W .)] ' (35) 
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Figure 2: Arrangement of wave vectors of all the modes resulted when the incident wave is of 
mode 2, ~£* , and when the anisotropy vector a has the direction as shown here. The grazing 
angle of the reflected mode 2, £ 2 , is less than specular one (specular direction is shown by 
broken arrow), and the grazing angle ipi of the reflected mode 1, £ ±, is even lower. The 
grazing angle tp of the transmitted wave ^ is even lower than tpi. We can imagine than at 
some critical value <p = (p cl the angle ipo becomes zero. It means that at ip < ip c i transmitted 
wave becomes evanescent and all the incident energy is totally reflected in the form of two 
modes. More over, there exists a second critical angle <^ c2 , when ipi = 0. It means that at 
p < (p C 2 the mode £ i also becomes evanescent. In this case all the incident energy is totally 
reflected nonspecularly in the form of mode 2. At the same time the two evanescent waves 
~£^o and £ 2 combine into a surface wave, propagating along the interface. The arrows over £ 
show direction of waves propagation with respect to the interface. In the figure there is also 
shown the basis which is used along the paper. It consists of unit normal vector n along normal 
(2-axis), unit tangential vector I (rr-axis) which together with n defines the incidence plane, 
and the vector t (y-axis) looking toward the reader, which is normal to the incidence plane. 



In the opposite case, when the incident mode is transverse one, the reflected mixed mode will 
have k 2 i± shown in (32). Therefore the difference A 2 i = (k 2 i± — &ir_i_)Ao\/^o is 



7/gsin(20 o ) + 2^(1 + ^)[1+V cos 2 (0 a )] -g 2 (l + 7/) 



v 2[i + ,;w(M ^-« 2 - < 3e > 

The changes of normal components with variation of 9 a according to (34), (35) and (36) for 
some values of dimensionless parameters 77 and q and vector a lying completely in the incidence 
plane, are shown in Fig. 1. From this figure it is seen that the strongest deviation of reflected 
wave from specular direction is observed for reflection of mixed to mixed mode. 

Since reflection of mode 2 is in general nonspecular, it can happen that the wave vectors of 
reflected and transmitted waves will be arranged as shown in fig. 2, and it follows that there are 
two critical angles for (p. The first critical angle (p c i (q 2 = 1/eo) is the angle of total reflection. 
At it the transmitted wave becomes evanescent. The totally reflected field contains two modes. 
At the second critical angle ip c2 , when q is in the range 

, (l + V )(l + rfcos 2 (9 a )) 
L<q < 1 + 77' ■ [6l) 

the reflected mode 1 also becomes evanescent. Together with evanescent transmitted wave 
the mode 1 constitutes a surface wave, propagating along the interface. In that case we have 
nonspecular total reflection of the mode £ 2 - 

In figure 3 it is shown how do the normal components of wave vectors change with increase 
of q, which is equivalent to decrease of tp. For eo = 1.6 the first critical angle corresponds to 
q ~ 0.8. The second critical angle corresponds to q — 1. 
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Figure 3: Dependence of dimensionless normal components of incident and reflected waves 
on q = kcostp /k ^/I^. The solid curve corresponds to the incident wave moving to the right 
kr(q) = k2r±/k ^/eo. The dotted curve corresponds to the reflected wave of mode 2 moving 
to the left kl(q) = k2i±/k 0y /e^. And the broken curve corresponds to the reflected wave of 
mode 1 moving to the left kl(q) = ku±/k ^/eo. It is seen that at q > 1 the mode 1 ceases to 
propagate. Its normal component kl(q) = y/1 — q 2 = —i\/q 2 — 1 becomes imaginary, therefore 
the reflected mode 1 becomes an evanescent wave. Together with transmitted wave, which 
becomes evanescent at q 2 = l/e , the mode 1 constitute the surface electromagnetic wave. 

3.2 Reflection and refraction from inside anisotropic medium 

The wave function in the full space is 

*(r) = 6(z < 0) ( e ^-^.+ £ Arrtp jl ^ il ) + e(z>0)e^ r (^ ej +^ mj ), (38) 

V j'=l,2 / 

where = £ + 'H, arrows show direction of waves propagation, ~$ • denotes the incident wave 

<— — >■ 

of mode j (j = 1,2), ip y (I = 1,2) denotes reflected wave of mode f, kj = (fey, kj r ±), 

k j> = (fey, —kjn±), k = (fey, ^Jkl — k 2 ), ipe^m, ~T^ e ,mj are fields and transmission amplitudes of 

TE- and TM-modes for incident j'-mode respectively. To find reflection ~j? and transmission 

r amplitudes (the arrow over them shows the direction of propagation of the incident wave 

toward the interface), we need to impose on (38) the following boundary conditions. 



3.3 General equations from boundary conditions 

Every incident wave field can be decomposed at the interface into TE- and TM-modes. In 
TE-mode electric field is perpendicular to the incidence plane, £ oc t, therefore contribution 
of j-th mode into TE-mode is (£j ■ t). In TM-mode magnetic field is perpendicular to the 
incidence plane, H oc t, therefore contribution of j-th mode into TM-mode is (Hj ■ t). For 
transmitted field in TE-mode we accept o e = t, l^, e = [kq x t], and for transmitted field in 
TM-mode we accept m = t, ~~£ m = —[k, x t]. 



3.3.1 TE-boundary conditions 

In TE-mode for incident j-mode we have the following three equations from boundary condi- 
tions: 

1. continuity of electric field 

(t ■ tj) + (t ■ t 1 )jf lj + (t ■ t 2 )J? 2j = 7> ej , (39) 



7 



2. continuity of magnetic field parallel to the interface 

(I ■ %) + (I ■ + (I ■ %)l?2 j = (I ■ [«o x t])-? ej = -Ko ± -? el , (40) 

3. and continuity of the normal component of magnetic induction, which for /j, — 1 looks 

(n • it,) + (n ■ + (n • ^ 2 )^ 2 j = (« • [«o x *])^ej = ko||T C j- (41) 

The last eq. (41) is, in fact, not needed, because it coincides with (39). It is a good exercise to 
check the identity of (39) and (41) using explicit expressions for £ and f£. Some examples of 
such substitutions will be presented in the Appendix A. 

3.3.2 TM-boundary conditions 

In TM-mode we have the equations 

1. continuity of magnetic field 

(t ■ %) + (t ■ + (* • %)^2j = (42) 

2. continuity of electric field parallel to the interface 

(I ■ + (l ■ + (l ■ t 2 )J? 2j = -(l ■ [k x t])-? mj = n 0± l> mv (43) 

3. and continuity of the normal component of field D 

(n ■ e£j) + (n ■ e£ 1 )'f lj + ( n ' s£ 2 )~^ 2 j = -( n ' i K o x A)~?mj = K o\\^mj- (44) 

Again we can neglect Eq. (44), because it coincides with (42), and again it is a good, though a 
little bit more difficult, exercise to check identity of these equations using explicit expressions 
for £ and H. In the following we will not show third equations like b23a and (44), because 
they are useless. 

Exclusion of !* e j from (39) and (40), and exclusion of l* m j from (42) and (43) give two 
equations for Pij, P 2j> which is convenient to represent in the matrix form 

'(* ■ %) + «o±(t • ((I ■ %) + «o ± (i • t 2 )) \ ( ~f x \ = _(il- %) + «o ± (i • 

K0±(t ■ %) ~ (I ■ tl)) («0±(t • %) ~ (I ■ t 2 )) ) U 2 J " V«0-L(* • %) ~ V • ~2j)J ' 

Solution of this equation is very simple if to take into account that reciprocal of an arbitrary 
2x2 matrix looks 

/ a b \ = 1 ( d -b\ , 4g s 
\c d J ad - be \ -c a J ' 

Therefore 

ftlj) = 1 ( ~ ( K 0±(t ■ %) - (I ■ t 2 j) ((Z • %) + K 0± (t ■ t 2 )) \ f(l ■ %) + K 0± (t ■ 

(47) 

Where D is determinant 

D = {{I ■ %) + Ko±(t • tl)) («o±(t • fe) - (I ■ t 2 )) - 

- ((I ■ %) + k ±(* • ^2)) («o±(* ' %) ~ (I ■ ^1)) • (48) 
Substitution of these expressions into (39) and (42) gives transmissions ~^ e ,mj 

GtJ = (|t • + ( t I* • fei ) (I2 ■) ' (49) 
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3.3.3 The most general case 

Above we considered the case when the incident wave has polarization vector ej with unit 
amplitude. (We remind that vectors ej are not normalized to unity.) To find later reflections 
from plain plates we will need a more general case, when the incident wave has both modes 
with amplitudes Xi^- To find amplitudes of reflected and transmitted waves in the general case 
it is convenient to represent the state of the incident wave in the form of 2 dimensional vector 

!*> = (£)■ <»> 

then the states of reflected and transmitted waves are also described by 2-dimensional vectors, 
which can be represented as 

|*>=(|)=it|*>, |^ )=(|»)=?V>, (5!) 
where TZ and T are two dimensional matrices 



We introduced the prime here and below to distinguish transmission and reflection from inside 
the medium from the similar matrices obtained for incident waves outside the medium. 

These formulas will be used later for calculation of reflection and transmission of plain 
parallel anisotropic plates. In the case of a plate we have two interfaces, therefore we need also 
reflection and transmission at the left interface from inside and outside the plate. Reflection and 
transmission from inside the plate can be easily found from symmetry considerations. Their 
representation is obtained from (47) — (49) by reverse of arrows and change of the sign before 
Ko±. After this action we find 

V P 21 V22 J V Vml Vm2 / 

Reflection from outside the medium is to be considered separately. 

3.3.4 Energy conservation 

It is always necessary to control correctness of the obtained formulas. One of the best controls 
is the test of energy conservation. One should always check whether the energy density flux 
of incident wave along the normal to interface is equal to the sum of energy density fluxes of 
reflected and transmitted waves, and the most important in such tests is the correct definition 
of the energy fluxes. In isotropic media it is possible to define energy flux along a vector n as 

/ j N (k-n) c e£ 2 + U 2 .... 

{J ' n) = — 7^^^' (54) 

(J-n^c^M. ( 55 ) 

47T 

In isotropic media both definitions are equivalent, because H = [k x £], and (k ■ £) = 0. 
The first definition looks even more preferable since the second one can be written even for 
stationary fields, where there are no energy flux. 

In anisotropic media only the second definition is valid, and because in mode 2 the field £ 
is not orthogonal to k, the direction of the energy density flux is determined not only by wave 
vector, but also by direction of the field £ itself. 
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3.4 Reflection and refraction from outside an anisotropic medium 



Let's consider the case, when the half space at z < is vacuum, and that at z > is an 
anisotropic medium. The incident wave falls onto interface from vacuum. The wave function 
in the full space can be represented as 



j'=e,m 



(56) 

where j,f denote e or m for TE- and TM-modes respectively, the term exp(i k\r)i() ■ with the 
wave vector A; = (fey, ^o_i_ = \Jkl — fcjj) describes the plain wave incident on the interface from 
vacuum. In TE-mode factor e = ~~£ e + contains e = t and = [7£ *]- m TM-mode 



factor t/> m = ~& m + T^m contains ^t e = t and e = — [ltot\. 



The reflected wave has the wave vector k o 



(fey, — fco±), and fields £ e = t, ^H e = [foot], 
— [Kot]. The refracted field contains two wave modes with wave vectors 



£m — t 7 and S m 

k i = (k\\,ki±), k 2 = (fc|| , A;2t-_l) an d electric fields ~£*i — e± — [alti] and £^ 2 = e 2 = 
a — it 2{alt 2)e(Q-^) / eq. Here k, = k/k, ki± = y^o^o — kp an d ^~2*± * s gi ven by (30). For 
incident TE-mode reflection p ee , p me and refraction Tj e amplitudes (j = 1, 2) are found from 
boundary conditions 

(ttj-tu + {tt 2 )^ 2e = 1 + y ee , (57) 

{l%)-t le + (Z^ 2 )^2e = "«0±(1 " T^ee), (58) 

{t%)-f le + {M 2 )^ 2e = 77 me , (59) 

(l^l)^le + (lt 2 )l>2e = -«0±^ me - (60) 

Exclusion of ee and ~~fF me leads to 



«0±(*^l) + (l^l)j f«0±(*^2) + (1^2) J J V^2 £/ 



2/t _L 

. . 



(61) 



and the solution 

m 



l 



(ko±(*^ 2 ) + (l^ 2 )) - («o±(*^ 2 ) - (i^ 2 ) 
- («o±(t^i) + (l^i)) (k ±(*^i) - (fill)) 



(62) 



where D e is determinant 

D e = (K 0± (tt 1 ) - (I^x)) (« ±(*^ 2 ) + (1^2)) " 
" («0±(*^2) " («£ 2 )) («0±(*^l) + (^1 

Substitution of "r^ e into (57) and (59) gives 



(t^0 (tt 2 ) \ (i> le \ _ ( i N 



(«t 2 ) 



In the case of incident TM-mode we have boundary conditions 

(^0^i m + {t!t 2 )-i> 2m = i + -jf mm , 

2m — «0J-U nun 



(63) 
(64) 

(65) 
(66) 
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(t~£ 1 )l> lm + (tt 2 )^ 2m = -f em , (67) 

(Ml^lm + {l%)~?2m = K ±J? em . (68) 



Exclusion of ~~jf me and p mm leads to 



V(*^i) + (^i)i r«o±(*^ 2 ) + (i^ 2 )] \ /^i m \ = /2« ±\ 



(69) 



Therefore 




(70) 



(l>lm\ = JL ( ( K *dtt 2 ) ~ {1%)) ~ («0±(t^ 2 ) + (^ 
V^ 2 J An ^ - (« ±(t^l) " (l^l)) (« ±(t^l) + (itl)) 

where D m = —D e (63). Substitution of Tj m into (65) and (67) gives 

(tern) _ ( (t^l) (t^ 2 ) \ ^ lB \ _ AA 

U m J"U^i) («t 2 ) J 1^2 J W' ( j 

In the general case, when the incident wave has an amplitude in TE-mode and amplitude 
if ' m in TM-mode, the state of the incident wave can be described by two-dimensional vector 

'MB- (72) 



and the states of reflected and transmitted waves can be represented as 



(73) 



where 7?. and 7" are the two dimensional matrices 



^=f I" V V (74) 

V Pme Pmm ) \ ^2e T 2m J 

4 Reflection and transmission of a plain parallel plate of 
thickness L 

Now, when we understand what happens at interfaces, we can construct [8] expressions for 
reflection, R(L), and transmission, T(L), matrices for a whole anisotropic plain parallel plate 
of some thickness L, when the state of the incident wave is described by a general vector lf )- 
To do that let's denote the state of field of the modes e.\ and e 2 incident from inside the plate 
onto the second interface at z = L by unknown 2-dimensional vector \x) (50). If we were able 
to find \x) we could immediately write the state of transmitted field 

T(L)|to> = rV>, (75) 
and the state of the field, reflected from the whole plate 

_a tr'tr 



R(L)\t ) = n\£ ) + T E(L)n | a»), (76) 
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Figure 4: Dependence of reflectivities |-R ee | 2 and |-R me | 2 of an anisotropic plate with eo = 1.6, 
i] = 0.8 and dimensionless thickness Lu/c = 10 on angle of the plate rotation around its 
normal, when the anisotropy vector a is parallel to interfaces and at <fi = is directed along 
fey. The incidence angle 6 is given by sin# = 0.9. 



where E(L), E(L) denote diagonal matrices 



E(L) = 



exp(iki±L) 






exp(ik 2i± L) 



E(L) = 



exp(iki±L) 






exp(ik 2r ±L) 



(77) 



which describe propagation of two modes between two interfaces. Here k±± = ^e G k { 
k 2r ± and k 2 i± are calculated according to (29) or (30) and (32), respectively. 
It is very easy to put down a self consistent equation for determination of 

|^> = E(L)T\£ ) + E(L)nE(L)n\lt). 



k?,, while 



(78) 



The first term at the right hand side describes the incident state transmitted through the first 
interface and propagated up to the second one. The second term describes contribution to the 
state of the | af) itself. After reflection from the second interface this state propagates to 
the left up to the first interface, and after reflection from it propagates back to the point z = L. 
Two terms at the right hand side of (76) add together, which results to some new state. But 
we denoted it and it explains derivation of the equation (76). 
From (76) we can directly find 



1^) 



E{L)KE{L)K] E(L)T\£o), 



and substitution into (75) and (76) gives 



T(L) 



T T 

- 1 ee - 1 em 

T T 

- 1 me -*■ mm 



T I- E{L)KE{L)K E(L)T, 



(79) 



(80) 



i(L) = 



Ree R 



1Z+TE(L)1Z I - E(L)H E(L)n 



E(L)T. 



(81) 



With these formulas we can easily calculate all the reflectivities and transmissivities for 
arbitrary parameters, arbitrary incidence angles, arbitrary incident polarizations and arbitrary 
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Figure 5: Dependence of transmissivities \T ee \ 2 and |T me | 2 of an anisotropic plate on angle <fi of 
the plate rotation around its normal, with all the parameters the same as shown in caption of 
fig.4. 

direction of the anisotropy vector a. In fig.4 we present, for example, reflectivities of TE-mode 
wave from a plate of thickness L such, that Lu/c = 10. The anisotropy vector is parallel to 
interfaces. Therefore, its orientation with respect to wave vector fc of the incident wave varies 
with rotation of the plate by an angle around its normal. The transmissivities of the same 
plate in dependence on the angle are presented in fig. 5. 

It is important to note that after transmission through anisotropic plate the transmitted ray 
has the same angle with the plate normal as the incident one, and it is interesting to investigate 
how the image of a real object splits, when observed through a birefringent ideal plain plate. 

5 Conclusion 

We have shown how to find polarizations of plain waves propagating in a single axis anisotropic 
media, how to calculate reflection and refraction at an interface between anisotropic and 
isotropic media, and how to calculate reflection and transmission of plain parallel transpar- 
ent anisotropic layers. We considered media without absorption, but inclusion of absorption 
does not provide any problem. 

We have shown that reflection from interfaces in anisotropic media is accompanied beam 
splitting, that reflection of the mixed mode is nonspecular and can be characterized by two 
critical angles. The first critical angle <p c i corresponds to total reflection with nonspecular 
double beam splitting. The second critical angle corresponds to total nonspecular reflection of 
mixed mode without beam splitting and creation of surface electromagnetic wave. The effect 
can be observed with the help of birefringent cone, when the incident beam is transmitted 
through the side surface and after reflection from the base plane goes out of the cone. 

We considered only one anisotropy vector, and it will be a good exercise to consider whether 
addition of a second anisotropy axis will not bring some new effects. 
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A Substitution of mode polarizations in boundary con- 
ditions for incident wave of mode 2 

A.l TE-mode 

Substitution of (15) for £ 2 , [a x n] for Z\ and (25) for 'H^ into (39) and (40) gives 

(t-a) + (t-[ax fci])-^ 12 + (t ■ a)^ 22 = l* e2 , (39a) 
k 2r± (t ■ a) - k 1± (t ■ [a x fei])"j^ 12 - k 2l± (t ■ a)~f 22 = k 0± !* e2 , (40a) 



where = (k\\,—ki±)/ki, k\ = k 0y /eo, k\± = ^Jk\ — /cjj, and k 2r j± are given in (30), (32). 
In equation (39a) we took into account that vector t is orthogonal to n, and in (40a) we used 
representation: n — K\\l + K±n, and the relation [I x n] = —t, which for any vector A gives 

(I ■ [k x A]) = ([I x k] • A) — - K± (t ■ A). (82) 

After exclusion of r e2 from these equations we obtain the equation 

(t ■ [a x 1^i})(k 1± + k 0± )~^ 12 + (t ■ a)(k 0± + k 2l± )^ 22 = (t ■ a)(k 2r± - k ±). (83) 

In the case of isotropic medium we have k 2r ± = k 2 i± = k±± = k±, — m — (fey, —k±)/k, 
therefore (83) is reduced to 

~ l" X ? r ^ tl2 + ^22 = PeO, ^e2 = (t ' a)r e0 , (84) 

(t-a) 

where 

(k± - k 0± ) 2k ± 

P e0 = Ti n — V' r e0 = 1 + p e0 = — — r (85) 

[k± + k ±) [k ± + fc _L) 

are the standard reflection and transmission amplitudes [5] of a pure TE-mode at an interface 
between isotropic media. 

If a = t we have the typical incident TE-mode, and p\ 2 is excluded because (£• [n r xo]) = 0. 

If (a • t) — we have the typical incident TM-mode, and from (39a), (40a) it follows that 

Pl2 = 

A. 2 TM-mode 

The similar substitutions of Tii^ into (42) for TM-mode gives 

^(t ■ [it 2 x a]) + ^(t • a)-? 12 + T^(* • [^2 x a])^ 22 = (42a) 

Kq Kq Kq 

For substitution of S into (43) we represent (15) in the form 

e 2 = [ K x[ax K ]] e ^-a A ^ ) , (86) 

where Ae(0) = ei(0) — e , and with account of (82) we obtain 

/. r -> n e i(#2r) Ae(# 2r ) v 
/t 2r _L(t • [K 2 x aj) (Z • a) - K 1± (t ■ a)p 12 - 

^0 e 
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- f« 2 ,±(* • [^2 x a]) 6 ^ + Ae M(Z • a)) y 22 = Kox^m- (43a) 

Exclusion of 7^ m with the help of (42b) gives 

{eok ±k 2r + k 1± h)(t ■ a)~ft 12 + ((e A;o±^2r + k 2 u_k 2 i){t ■ [1c 2 x a]) + k 2 Ae(0 2 i)(l ■ a))~^ 22 = 

= - e k ±)(t ■ [it 2 x a]) - £; 2 Ae(# 2 r)(J • a). (87) 

Together with (83) we have two equations for determination of p 22 and p 12 . 



In the case of isotropic medium we have k 2r ± = k 2 i± = k\± = k±, k 2r = k 2 \ = ku = k, 

Ae = 0, e-i — e , lt 2 = it = (k\\,k±)/k, 1c 2 = 1c i = 1c = (k\\,—k±)/k. Therefore (87) is 
reduced to 



(88) 



(t ■ a)-f 12 + (t ■ [1c x a])^ 22 = (t • [it x a]) ^ e °f 0± | = (t • ["^ x a])p m0 , 

+ e fco±J 

and solution of (88) with (84) gives 

v (t ■ a) ((t ■ [it x a])p m0 - (t ■ [fe x q])p e0 ) 

//l2 " (t • [fe x a]) 2 + (t • a) 2 ' m 

v (t ■ [1c x q])(t ■ g x a])p m0 + (t ■ a) 2 Peo (Qn , 
(f22 ~' (t ■ [1c x a]) 2 + (t • a) 2 • (9Uj 

For r m it follows from (42a) that 

2k 

■? m = (t-[ltx a])(l + p m0 ) = (t • [7c" x a])- ■ — = (t ■ [it x a])r m0 , (436) 

fc_L + e k ± 

where p m o and r m o are the standard reflection and transmission amplitudes [5] of a pure TM- 
mode at an interface between isotropic media. 

Thus, for an internal reflection from an interface of isotropic dielectric with vacuum we 
obtained reflection and refraction amplitudes for an incident field with an arbitrary polarization 
e = [k x [a x k]], which is determined by some unit vector a. If a = t, then 7^12 = an d 
7^22 — Peo- If ( a 't) — 0, then again ~~/? 12 = 0, but 7^22 — Pmo(t- [it x a])/(t- [1c x a]). If a = 1c , 
then 7^22 1S divergent, however the reflected field £ 2 ~$ ' 22 has the finite value p m o(t • [it x a]). 



B History of submission and rejection 

The paper was submitted to Am.J.Phys on September 1 of 2010. It was rejected on September 
24 because of negative reports of two referees. The first referee said that he is lazy to read the 
manuscript with pencil, but he saw that sections 2 and 4 are absolutely not needed, because 
everything about reflections is much better said in textbooks by Jackson and Griffith. The 
second referee rejected because we, he said, erroneously told that the recent paper was publish 
long ago at 1977. He said that since then there were many papers on optical reflection and 
transmission. 

If we could reply to referee we would mention that in textbooks by Jackson and Griffith 
there are no word on anisotropic media. The referee overlooked the main point of our paper. 
As for claim of the second referee, we would like to say, that he can try to seek on AJP home 
page a paper with key words "electromagnetic waves in anisotropic media." Then he will find 
the first article published in 1977. So the second referee also overlooked the main point of our 
article. 
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